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D ' Abstract. We derive a generalized linear dispersion relation of waves in a strongly 

y . magnetized, compressible, homogeneous and isotropic quasineutral plasma. Starting 

O ' from a tvi^o fluid model, describing distinguishable electron and ion fluids, we obtain 

c/3 , a six order linear dispersion relation of magnetized waves that contains effects due to 

C/2 ' electron and ion inertia, finite plasma beta and angular dependence of phase speed. 
We investigate propagation characteristics of these magnetized waves in a regime 

^ ' where scale lengths are comparable with electron and ion inertial length scales. 



This regime corresponds essentially to the solar wind plasma where length scales, 
Q,' comparable with ion cyclotron frequency, lead to dispersive effects. These scales in 

conjunction with linear waves present a great deal of challenges in understanding 
the high frequency, small scale dynamics of turbulent fluctuations in the solar wind 
Kj^ . plasma. 

(N 

O 

J-/ , 1. Introduction 

^^ ' Solar wind plasma is an admixture of waves, structures, and turbulent fluctuations 

^^ . that comprises multitude of length and time scales. Owing to a great deal of dis- 

parity in the time and spatial scales, solar wind plasma exhibits rich and complex 
dynamical evolution. Since the solar wind plasma is strongly magnetized, the pres- 
k> ■ ence of waves and their consequent interactions with fluctuations complicate our 

'V^ . understanding of many aspects. For instance, the solar wind fluctuations yield a 

C^ ' composite spectrum [1] as described in the schematic of figure 1. This spectrum 

describes the power spectral density (PSD) as a function of frequency and can be 
divided into five distinct regions. The frequencies smaller than 10^ Hz, namely re- 
gion I, lead to a PSD that has a spectral slope of 1 [1, 2]. The region II extends 
from 10^ Hz to or less than ion/proton gyrofrequency where the spectral slope 
exhibits an index of 3/2 or 5/3. The latter, a somewhat controversial issue, is char- 
acterized essentially by fully developed turbulence. The region II connects to the 
region HI by means of a spectral break at lengthscales corresponding to ion inertial 
lengthscales and frequencies corresponding to ion gyrofrequencies. The onset of the 
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spectral break is disputed. This regime is often referred to as dissipative regime 
which exhibits a PSD with a much broader spectral slope that varies between 2 
and 5 [6, 3, 4]. Notably, the dynamics of the lengthscales in this region cannot 
be described by the usual compressible or incompressible MHD models that pos- 
sess characteristic frequencies smaller than ion gyro frequencies. A two fluid MHD 
model needs to be invoked to examine the dynamics of these high frequency and 
fast time scale processes. What is notable in these spectra is the waves in differ- 
ent regimes that interact with the turbulent fluctuations and influence the cascade 
dynamics [f6, f7, f8]. Despite their prominent existence in different regimes, less 
is understood about the evolution of these waves and their dynamical role in gov- 
erning the turbulent spectrum. For example, it is argued that the spectral break 
results from the damping of ion cyclotron waves. The latter is contrasted by Shaikh 
& Zank [5] who report that Hall forces could play a critical role in producing the 
spectral break. Similarly, the excitation and interaction of whistler waves in high 
frequency turbulence has been disputed recently [7, 8, 9, fO]. Motivated by these 
issues, we in this paper investigate linear waves in strongly magnetized plasma that 
are relevant to the understanding of the solar wind turbulence spectra depicted in 
Fig f. Our main objective here is to develop a comprehensive understanding of the 
propagating linear waves in strongly magnetized plasma especially in the regimes 
HI, IV and V of Fig I. We further concentrate on a high beta plasma, because the 
solar wind plasma fluctuations are characterized typically by a high beta plasma 
beyond 1 AU (astronomical unit). 

In section 2, we describe our two fluid model that contains all possible modes 
in the magnetized solar wind plasma. A generalized linear dispersion relation is 
derived. In section 3, we describe various possible roots, corresponding to modes. 
This section also describes the solution of the generalized dispersion relation and 
focuses on waves especially associated with a high beta (pressure and magnetic en- 
ergy ratio) plasma. Section 4 deals with the effect of plasma beta on the propagation 
characteristics of dispersive waves. Finally, a summary is presented in section 5. 

2. Linear Dispersion Relation for two fluid plasma with flnite 
electron mass 

We start with a two- fluid model (ion with suffix i and electron with suffix e). 
Using usual notations and with finite electron mass, i.e., TOg ^ 0, the ion and 
electron continuity equations, momentum equations, energy equations together with 
the Maxwell's equation (neglecting displacement current) are written as (in CGS 
Gaussian unit) 

^ + V.(n,uO-0, 11^ + V • (neUe) = 0, (2.1) 

miUi [ — -^ + u; • Vui ) = -Vpi + TiiC ( E + -Ui X B ) , (2.2) 



dt 



c 



rUeUe ( -^ + Ue • VUe j = -Vpe - UeC f E + -Ue X B ) , (2.3) 

-^+\li-ypi+ J,p,V • Ui == 0, ^ -F Ue • Vpe + 7ePeV ■ Ue = 0, (2.4) 

ot ot 
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Figure 1. Schematic of the solar wind composite spectrum as a function of frequency 
(wavenumber) . Magnetized waves play a critical role in the transition from region II (MHD 
regime) to region III (kinetic or Hall MHD regime). The onset of regions HI and IV is asso- 
ciated with high frequency whistler waves where electron and ion motions are decoupled. 
Regions IV and V are identified as whistler wave cascade regimes. 



Aire 

V X B = (n^Ui - rigUe), 

c 

5B 

V-B = 0, 



(2.5) 

(2.6) 
(2.7) 



where UijU,, are the ion and electron number densities, mj,me are the ion and 
electron mass, u,, Ue are the ion and electron velocity, Pi,Pe are the ion and electron 
pressure, 7i,7e are ratio of the specific heats for ions and electrons, respectively, and 
E is the electric field and B is the magnetic field. We use the linearized forms 
of the above equations with the following perturbation scheme for any variable 
ipa = i'ao + i^a, (a = i,e): 



naO + fia, Pa^PaO+Pa, 



B — Bq + B E — E; 



(2.8) 



with riaOjPaOT Bq are constant and uniform in space. We assume the space-time 
dependence of a perturbed variable ip{x, y, z, t) as, ip ~ exp[i(k-r — wt)]. Further, we 
introduce the plasma displacement vector ^^ as, Uq = d^a/dt — —iu)S,ai with ^^ as a 
vector. From the linearized equations for ion and electron continuity equations, and 
the energy equations, one gets, (replacing d/dt — > — icj; V — > ik) the following: 



-inao(k • ^a), Pa = -i"taPao{^ ' ^a), 



(2.9) 
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The linearized form of (2.5) after combining with Unearized form of (2.6), can be 
written as, (with riio = Ueo = no), 

kxkxE=-c.2l^(^i-^,) (2.10) 

Finally, we write the linearized forms of momentum conservation equations as: 



lUJ 



uj m^noCi == -"ftPtokik ■ ^i) + uqc E ^i x Bq , (2.11 



c 



ILU 



-u mena^i = -7ePeok(k • ^e) + noed,, k x k x E - uqc f E ^e x Bq 

(2.12) 

where we used Uio = Ueo = uq and replaced ^e the displacement vector for electron 
by the relation obtained from (2.10), i.e., 

Ce = Ci + -, J k X k X E; 

and introduced the electron incrtial length de ~ c/ujpe, bjp,, being the electron 
plasma frequency. From eqn. (2.10)taking a dot-product with k, and also using the 
continuity equation (2.9) we get 

k-$i = k-^e, ne = ni, (2-13) 

where the second equation denotes the quasineutrality condition. Equations (2.11) 
and (2.12) are our starting equations for deriving the dispersion relation for a two 
fluid plasma. We followed the procedures adopted by Ishida et al., to get the dis- 
persion relation and we omit the derivation in this work for brevity. We proceed by 
adding both the equations, eliminating the perturbed electric field E in terms of ^i 
using (2.11) and finally taking the dot product of the resulting equation for ^i with 
three independent vectors, k, b and k x b, where b is the unit vector along the 
unperturbed magnetic field Bq. The final generalized dispersion relation reads as 

[cj2(l + k^dl) - fcVj cos^ 0] [cj4 - fc2(vj + Cl)uj^ + k^Vlk^Cl cos^ 6 
+uj^k^dl{uj^ - k^Cl)] ^ up-Vlk^d\k^ cos^ Qiu? - k^Cl); (2.14) 

where di — c/upi, ujpi being the ion plasma frequency, and is the angle between 
the wave vector k and unperturbed magnetic field Bq, and Va, Cs are the Alfven 
speed and sound speed, respectively. This dispersion relation (2.14) incorporates 
both the ion and electron inertial length scales, di and de, respectively. Moreover, 
equation (2.14) is essentially the same as the dispersion relations obtained by earlier 
authors in Refs [13, 15, 14]. It is easy to see that for electron mass irie — ?■ 0, which 
is equivalent to putting d^ = 0, eqn. (2.14) reduces to 

y - k^vi cos^ e] [lj^ - k^ivi + cl)iu^ + k^vlk'^cl cos^ e] 

= ijj'^Vlk^djk^ cos^ 6'(cj2 - k^Cl); (2.15) 

which is the usual Hall MHD dispersion relation obtained by Ishida et al. [12]. For 
both di =Q, de = eqn. (2.14) reduces to the dispersion relation for an ideal MHD, 
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showing the existence of shear Alfvcn wave, fast and slow magnetosonic waves. 

From the dispersion relation (2.15) Ishida et al [12] have shown some interesting 
effects of the presence of ion inertial length. It is shown that for the limit kdi > 1, 
the incompressible MHD Alfven wave becomes compressible and the MHD com- 
pressible slow wave becomes incompressible. While we postpone similar studies for 
the inclusion of electron inertial length for our future work, we shall study the effects 
of plasma (3 on the linear waves in solar wind from the dispersion relation (2.14). 
For that we normalize the frequency cv by some characteristic frequency cjo,so that 
the dimensionless frequency w = uj/ujo = w. Let Uq be some characteristic speed. 
Thus Lq — Uq/ujo will be some characteristic length. Choosing ujq = (/cVa), and 
note that 



C| ^ IPo/Po ^ 47r7Po Snpa ^ Ml 

VI Bl/Airp, Bl "" El ^' Ml 



where Ma and Ms are Alfvenic and sonic Mach numbers. 

We next normalize w by kVA and Lo/k ~ Vp, the phase velocity, and uj/kVA is the 
normalized phase velocity. The linear phase velocity relation can then be expressed 
as 

[^2(1 + (kdef) - cos^ e] = [v^ - (1 + p)vl + (3 cos^ 6 + vl{kd,f{vl - /3)] 

= vlikd.Y cos^ 6{vl - 13) (2.16) 

So the final form of the dispersion relation given by eqn.(2.16) is essentially an 
equation showing the angular dependence of normalized phase velocity to plasma 
f3 and the ion electron inertial length scales. 

3. Dynamics of linear dispersive waves 

We have developed a matlab code to solve the generalized dispersion relation Eq. 
(2.16). The dispersion relation is 6th order in frequency. Hence six roots are ex- 
pected. Three roots correspond to forward and the remaining three represent back- 
ward propagating waves. The generalized dispersion relation, Eq. (2.16), spans a 
wider parameter regime and exhibits a variety of waves for various angle, beta, 
length-scales comparable with the electron and ion inertial skin depths. In this pa- 
per, we nonetheless restrict ourselves to the parameter regime that is a representive 
of the solar wind plasma corresponding to the length scales that are associated with 
the ion cyclotron frequency. 

As a first step, we verify the consistency of our equation by comparing it with the 
MHD waves 11. For this purpose, we use dt = de = that reduces the dispersion 
relation to the usual compressible MHD relation. We retrieve MHD waves that are 
consistent with Ref 11. This is shown in Fig (la) that describes a phase velocity 
variation of the MHD waves as a function of angle between the propagation wave 
vector and the mean magnetic field, i.e. k\\ = kcos9 where fcy = k ■ Bq. We find 
three forward propagating modes with positive phase velocity that co-exist with 
three backward propagating modes with negative phase velocity. It is noteworthy 
that the shear Alfven and slow modes are partially overlapped. In this regime, linear 
dynamics is entirely governed by the magnetosonic waves which is shown by the 
top curve in Fig (la). A pure magnetic perturbation propagating orthogonal to the 
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(a) MHD regime 
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Figure 2. (a) Linear waves in MHD regime. Forward propagating fast, slow, and Alfvenic 
modes are shown. Tlie top curve describes Alfvenic mode, whereas the remaining two 
curves are overlapped and represent fast and slow modes. Our results are consistent with 
Ref [11]. (b) High frequency waves in kdi > 1 and kde < 1 (c) kdi > 1 and fcde > 1 
regimes. It is clear from (b) and (c) that the propagation characteristic is modified dra- 
matically in the short scale regimes. The phase speed increase for shorter scales. The 
parallel propagation is faster than the oblique one. 



constant magnetic field in this regime behaves electrostatically and tends to move 
as a magnetosonic mode. 

We next investigate the effect of a finite plasma beta on the propagation of the 
small scale magnetized waves to explore their dynamics in regimes III and IV in Fig 
(1). We choose parameters (3 = 3,k — 0.5, de = 0.1, di — 43de that are characterize 
the high frequency modes between regions II and III. Since this parameter regime is 
not drastically different than the typical MHD regime, a non trivial modification in 
the dispersion of MHD waves is expected. For instance, the linear dynamics in this 
regime is governed predominantly by the higher frequency and short scale waves, 
namely kinetic Alfven and magnetosonic waves. The density perturbation becomes 
significant in this regime. The dispersion of high frequency kinetic Alfven, shown in 
Fig (2b), differs significantly from that in Fig (2a). The phase velocity of the kinetic 
Alfven waves is increased, while the obliquely propagating kinetic Alfven modes 
continue to remain unaffected. The shear Alfvcn/slow modes decouple clearly. The 
fast modes travel with higher phase speed, whereas the slow modes propagate slower 
(than in Fig 2a). It is further clear from Fig (2b) that the propagation of the 
fast/slow modes depends crucially on the alignment of their wave vector relative 
the mean magnetic field. For a highly oblique propagation, these waves hardly move. 
Figure (2c) describes a regime where kdi > 1 and kdg ^ 1. This is a regime where 
MHD modes are drastically altered and low frequency whistler modes start to play 
critical role. Shown in Fig (2c), the top curve is the whistler branch that survives 
whereas the bottom curves are reminiscent of MHD modes (Alfven and fast/slow) . 
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(a) kd,=4.5,kd^=0.1 (b) kd,=45,kd_=1 .0 

20 





Figure 3. Figures (a) and (b) are shown for a large k. Note that the phase velocity of 
the small scale wave is consistent with Fig (2). 



(a) Low beta : kcl.=4.5,kd =0.1 




(b) High beta : kd.=4.5,kd =0.1 




Figure 4. Effect of plasma beta is shown for large scale waves, (a) describes a low 
plasma beta regime, whereas (b) represents the higher beta regime. 



As seen in Fig (2c), the whistler modes have higher frequency and phase speeds. 
They propagate predominantly along the field lines, whereas oblique whistlers have 
smaller phase speeds. At an angle 6 — 7r/2, the whistler modes do not propagate. 
They are transformed into high frequency electrostatic modes. Figure (2c) describes 
waves that possess length scales smaller than the ion inertial length, but bigger than 
electron inertial length scale. We can retrieve the small scale waves by choosing 
the small characteristic length scales of these waves relative to the ion inertial 
and electron inertial length scales. This is shown in Fig (3). This figure is further 
consistent with Fig (2) where the dynamics of small scale is described by choosing 
a relatively small magnitude of k compared to di and de lengths. 
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(a) Low beta : kd-45,kd =1 .0 



(b) High beta : kd-45,kd =1 .0 
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Figure 5. Phase speed variation with respect plasma beta for small scale dispersive 
waves, (a) and (b) describe respectively low and high beta regimes. 

4. Effect of plasma beta on the dispersive waves 

We next investigate effect of plasma beta (ratio of plasma pressure and magnetic 
energy) on the propagation of dispersive waves in the solar wind plasma. We keep 
the angle of propagation constant (9 ~ 45°) and vary plasma f3 and study the waves 
corresponding to the small and large scale modes. A snap shot for the large scale 
waves is shown in Fig (4) that depicts effects of low (Fig 4a) and high (Fig 4b) 
plasma beta on the propagating waves. We find that the large scale Alfven and 
slow waves continue to remain unaffected with regard to the beta varation. The 
propagation speed of the fast waves appears to grow linearly with the plasma beta. 
Hence there is a dramatic speed enhancement occurs for the fast waves. This is 
shown in Fig (4b). 

Interestingly, we find that the propagation property of the dispersive waves cor- 
responding to the small scale modes does not depend critically on the plasma beta. 
Almost all of these small scale waves remain unaffected relative to the beta varia- 
tion. This is shown in Fig (5a) for the low beta regime, whereas Fig (5b) describes 
the effect of high plasma beta on the small scale dispersive waves. 



5. Summary 

A major outcome of work is the inclusion of electron and ion inertia in strongly 
magnetized waves associated with a high beta plasma. We invoked a two fluid or 
Hall MHD model to examine the dynamics of these high frequency and fast time 
scale processes. We see that for the length scales kde ^ 1, where there the turbulent 
spectra shows a different power law behavior, the linear waves are significantly 
modified. It is clear from our work that the evolution of the linear waves in the 
MHD regime is modifed significantly in the dissipative or dispersive regime where 
kde > kdi > 1. This is the regime where solar wind turbulence is affected non 
trivially by the linear collision less waves. We find from our linear theory that MHD 
modes, such as Alfven and fast/slow waves, are inconsequential in the small scale 
kde > kdi > 1 regime. This is because the MHD modes are relatively large scale and 
low frequency modes in comparison with the kinetic Alfven and/or whistler modes 
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that are typically excited in the kde > kdi > 1 regimes of Fig (1). Because of the 
great disparity in the length and time scales, the waves in the two distinct regimes 
(i.e. regimes III, IV and V in Fig 1) do not couple efficiently. Consequently, the linear 
dynamics of modes in the regimes IV and V is governed predominantly by fast and 
whistler modes respectively. The latter have a greater phase velocity relative to the 
MHD modes. Our results are important particularly to understand linear waves 
in the solar wind turbulence in the regime where small scale fluctuations exhibit 
dispersive characteristics [1, 2, 4, 6, 5, 8, 9, 10]. 
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